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Abstract 

This paper discusses multi-Skyrmions on the 3-sphere S s with variable radius L using 
the rational map ansatz. For baryon number B = 3, . . . , 9 this ansatz produces the 
lowest energy solutions known so far. By considering the geometry of the model we 
find an approximate analytic formula for the shape function. This provides an insight 
why Skyrmions have a shell-like structure. 



PACS-number: 12.39.Dc 



27| pages, 14 figures 



* email: S.Krusch@damtp.cam.ac.uk 



1 Introduction 



In this paper we discuss new developments in the SU(2) Skyrme model |lf and its gen- 
eralization to the case where the physical space is a 3-sphere of radius L. The Skyrme 
model is a nonlinear field theory of mesons whose field configurations are labelled by an 
integer, the topological charge. This charge can be identified with the baryon number 
B [|]]]. Static field configurations which minimize the energy of the Skyrme model for a 
given baryon number B are called Skyrmions. When the theory is quantized, the Skyrme 
model not only describes the proton and the neutron reasonably well |2|, |3| but also larger 
nuclei 0, [|. However, in order to be able to perform the quantization it is important to 
reach a good understanding of the classical solutions. In Skyrmions were calculated 
numerically for small baryon number B, and it was shown that they have certain discrete 
symmetries. These symmetries have been confirmed by the rational map ansatz p|] which 
also reproduces the energies of the Skyrmions with good accuracy. 

From a mathematical point of view, field configurations in the Skyrme model are rep- 
resented by maps from (IR 3 U {oo}) = 5 3 to 577(2) = 5 3 . Therefore, it is natural to 
generalize the model such that physical space is a 3-sphere of radius L. From a physical 
point of view, a Skyrmion on a 3-sphere describes a finite baryon density H. Reducing 
L increases the baryon density. Varying L the 5 3 model exhibits phenomena such as 
localization-delocalization transitions and the restoration of chiral symmetry. These re- 
sults can be compared with Skyrme crystal calculations ||, [K| |XXf| , but are also interesting 
in their own right. 

In this paper we only consider static configurations and their energies. In particular, we 
generalize the rational map ansatz to describe Skyrmions on 5 3 . The energies of our ansatz 
are the lowest ones known so far. Moreover, these approximations have a well defined limit 
for L — > oo, namely, the rational map Skyrmions in flat space. Geometric considerations 
lead to an analytic ansatz for the shape function which is completely specified by a small 
set of parameters. We show that this ansatz agrees well with the numerical solution and 
captures the behaviour of Skyrmions on 5 3 . 

In the following section we review the Skyrme model on general 3-dimensional mani- 
folds. The first part focuses on the geometrical meaning of the Skyrme energy. The next 
part clarifies the relationship between the geometric formulation in flat space and the stan- 
dard formulation in terms of the Lie group SU(2). In section ||] the model is generalized 
to the 3-sphere. We describe the rational map ansatz in detail and also recall the doubly 
axially-symmetric ansatz | I2| . In section [| we derive an analytic ansatz for the shape 



function and discuss its symmetries. For special values of B the energy can be calculated 
explicitly as a function of the radius L. In section [| we calculate the shape function numer- 
ically. We also compare the numerical shape function to the analytic shape function of the 
previous section. Finally, we use the analytic shape function to approximate Skyrmions in 
IR 3 . In section 13 we discuss chiral symmetry and phase transitions. 
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2 Geometry of Skyrmions 



In this section we describe static solutions of the Skyrme model on general 3-dimensional 
manifolds. We follow Manton's approach || and present a geometric point of view. First, 
we introduce the geometric notion of the strain tensor and construct the Skyrme energy for 
general manifolds. Then we discuss the properties of the energy density and derive some 
formulae which will be important in the following sections. Finally, we show that in flat 
space the geometric energy density is equivalent to the standard energy density. 

A field configuration is a map tc from a physical space S to a target space E. Both S 
and E are 3-dimensional Riemannian manifolds which are connected and orientable, and 
their metrics are and T a p, respectively. We denote by x a point in S and tt(x) its image 
in £. Choosing dreibeins e m l (x) on S and ^ a (n(x)) on E, we can define the Jacobian 
matrix 

J m ^x) = e m l (x)(d t 7r a (x))(^(n(x)). (2.1) 

The matrix J mM (x) is a measure of the deformation induced by the map ir at the point x 
in S. However, it is not unique since a rotation of the dreibeins (^(-ir^x)) does not affect 
the geometrical deformation. This leads us to define the strain tensor 

D mn = (JJ T ) m n = e^id^e^id^)^ (2.2) 

which only depends on the metric on E. Here and in the following we suppress the direct 
reference to x. The strain tensor is symmetric and positive semi-definite[] but it is not 
invariant under rotation of the dreibeins e m \ In fact, under an orthogonal transformation 
O, D transforms into T DO. A well known result from linear algebra is that the char- 
acteristic polynomial -P(A) = det(D — XI) is invariant under orthogonal transformations. 
Denoting the non-negative eigenvalues of the strain tensor D by Af , X 2 ,, and A| we obtain 
the following invariants of D: 

Tr D = X 2 + X 2 + X 2 

|(Tr Df - |Tr D 2 = AfA| + A|A| + X 2 X 2 3 (2.3) 
detD = X 2 X 2 X 2 . 

As we will demonstrate in this section, choosing the energy of a field configuration to be the 
integral of the following sum of these invariants generalizes the Skyrme model to arbitrary 
manifolds S and E: 

E = [ Xl + X\ + Xl + X\X\ + X\X\ + XlXl (2.4) 

where the integration measure is V det t d 3 x. The first three terms are only quadratic in 
the derivatives and correspond to the "sigma model term". The configurations which are 

1 Generically, J is non-degenerate so that D is positive definite. However, there are submanifolds of 
zero baryon density, i.e. det(J(a;)) = 0. For further discussion see [O. 
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minimizing this term are known as harmonic maps |TJj]. Geometrically, n induces a linear 



map 7r* mapping the unit vector e m l which is formed by the mth dreibein to the vector 
e m J <9j7r a . The sum of the squared lengths of these vectors is the "sigma model term". Thus, 
only the lengths of the vectors 7r*(e m J ) are important. 

The geometric meaning of the quartic terms can be understood by considering the area 
element e qmn e m t e r J formed by two dreibeins. This area element is mapped to the following 
area element in S: 

e^ejejd^d^. (2.5) 

The sum of the squares of these area elements is proportional to the "Skyrme term" . This 
term is also a measure of the angular distortion of the map n*. 

So far only the first two invariants of the strain tensor have been used. The third 
invariant is the square of the determinant of the Jacobian, (det J) 2 . Since the manifolds 
S and E are orientable, det J is globally well defined so that we can set %/ det D = det J. 
det J locally changes the integration measure on S into the measure on S. Therefore, 
the integral of det J over S is equal to the volume Vol(E) times the degree of the map 
7r. The degree of a map is a topological invariant which, roughly speaking, measures how 
many preimages a point tt(x) of S has in S, and hence takes integer values. Any field 
configuration can be labelled by its degree. In the Skyrme model the degree of the map n 
is identified with the baryon number B, and we obtain 

We will call classical field configurations which minimize the energy (|2.4j) Skyrmions. It 
is worth pointing out that there is an alternative expression for equation ( |2.4| ). By "com- 
pleting the square" we can rearrange the Aj in the following way: 

E= f(X 1 ± A 2 A 3 ) 2 + (A 2 ± A3AO 2 + (A 3 ± AiA 2 ) 2 =F 6A!A 2 A 3 . (2.7) 

Either the upper or the lower signs are chosen such that the integral over the last term in 
equation ( |2.7| ) is non-negative. Applying equation ( |2.6| ), this integral is just 6 \B\ Vol(S). 
Therefore, the energy is bounded below by the so-called Faddeev-Bogomolny bound [|TJJ: 

E > 6 \B\ Vol(S). (2.8) 

Equation ( |2.7|) also gives rise to two sets of three Bogomolny equations which are satisfied 
if and only if the Faddeev-Bogomolny bound is attained, 

Ai = TA2A3, A 2 = TA3A1 and A 3 = ^AiA 2 . (2.9) 

The only non-trivial solutions of both sets of equations are Af = A 2 = A| = 1, and it follows 
that B = AiA 2 A3 = =pl. Therefore, the strain tensor is the identity map, and the map n is 
an isometry. For S = SU(2), this case can only occur if the physical space is a 3-sphere of 
radius L = 1. This has already been proven in [ITH]. 
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The relative importance of the "sigma model term" and the "Skyrme term" can be seen 
from their scaling behaviour. For this purpose we consider a family of metrics L 2 tij where 
L is a constant length scale. The dreibeins e m * are, informally speaking, the square roots of 
the inverse metric ft of S. Therefore, they are proportional to r so that the eigenvalues 
Af of the strain tensor Dij are proportional to jp. Since the measure ->/det(L 2 t^) d 3 x 
scales with L 3 , the "sigma model term" scales like L, whereas the Skyrme term scales like 

For large radius L one might expect the "sigma model term" to be dominant. Yet, this 
is not the case because, as we shall see, the Skyrmion becomes localized and — just as in 
flat space — both of these terms are equally important. However, for small radius L the 
configuration is delocalized, and the "Skyrme term" will become dominant. 

In the following, we relate this geometric formulation to the standard Skyrme model in 
flat space. For the remainder of this section S will correspond to the flat space M. 3 and E 
to the Lie group SU(2) = S 3 . The basic field is the SU(2) valued field 

£/(x) = cr(x) + mt(x) • t (2.10) 

where r« are the Pauli matrices which are Hermitian and satisfy the algebra TjTj = 6^ + 
itijkTk- Since U is an element of SU(2), the fields a and it obey the constraint cr 2 +7T 2 = 1. 
It is worth noting that tt in ( |2.10| ) play the same role as the coordinates Hi in the discussion 
above, and a is just a function of those coordinates. The static solutions of the Skyrme 
model can be derived by varying the following energy fTj: 

E = j r-^Tr(^-^([^,^][^,^])jd 3 x (2.11) 

where Ri = (diU)U^ is an su(2) valued current. A static solution of the variational equa- 
tions could also be a saddle point. Only solutions which minimize the energy are called 
Skyrmions. In order to show that the energies (|2.4j) and ( 2.11|) are equivalent, we first 
relate the strain tensor to the current Ri. Since a is determined by 7r we can calculate 
the induced metric r a/ g on the target space £ in terms of the fields tt: 

Ta/3 = Oa/3 + — jj- . (2.12) 

Starting from equation (|2.2j ) and noting that in flat space the dreibeins can be chosen to 
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be Kronecker deltas we obtain the following expression: 



D 



ij 



a/3 



(2.13) 



diodjO + di'K ■ dj7r 



(2.14) 



-Trma + iT-d^ma-ir-d^)) 



(2.15) 



(2.16) 



--Tr(i^) 



(2.17) 



Equation (|2.14|) follows from (|2.13|) by using the chain rule and the formula for the metric 
( |2.12j ). Conceptually, the step from equation ( j2.14| ) to equation ( j2.15| ) is very important. 
The metric r a p is expressed with the group multiplication and the trace which provides 
a scalar product in su(2). This is the transition from geometric to Lie group language. 
Equation ( |2.16|) is a trivial consequence of definition (|2.10|) , and the last equation follows 
from UW = 1. Now, we can show that the two energy densities are equal. Since Ri is an 
su(2) current, and therefore anti-Hermitian, it can be expressed as Ri = Rf ir a . We start 
with the energy density in ( p. 11 ) and expand Ri in terms of Pauli matrices: 



)-— Tr 
; 16 


[r^R^, ir v ]Rf '/#[iv, iv]) 


(2.18) 


RtRt + \ 


Tr ^ £ ^fif^6 W v) 


(2.19) 


= D U - 


+■ \ (flWi^ ~ RiRjRiRj) 


(2.20) 




-Da ~\~ — (yDaDjj DijD ji) 


(2.21) 




= Tr D + - (Tr D) 2 - -Ti D 2 . 

2 K ' 2 


(2.22) 



Equation (|2.19|) follows from (p.18 ) by inserting the commutation relations of the Pauli 
matrices T. In the following equation we use the contraction e^eum = $ji$km — Sj m Ski- In 
equation (|2.21|) we apply the result ( |2-17| ). The last equality follows from the definition of 



the trace. Equation fl2.22| ) is the energy ( |2.4j ) because of ( |2.3| ). Therefore, we have shown 
the equivalence of the two approaches for S = M. 3 . It is worth noting that in the geometric 
picture the "Skyrme term" is related to the square of an area element, whereas in the 
standard approach it depends on the structure constants of the Lie algebra su(2). Both 
interpretations lead to generalizations. While in our geometric interpretation it is natural 
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to consider different 3-dimensional manifolds for S, and maybe for E, the Lie algebra 
approach leads to generalizing the Lie group SU(2) for example to SU(N). 



3 Skyrmions on S* 3 and Rational Maps 



In the previous section we have established the equivalence of the two energies (|2.4j ) and 
( P- 11| ) m IR 3 - In this section we will use ( j2.4| ) to generalize the model to S = such that 
from now on physical space is a 3-sphere of radius L. This includes the original model if 
we take the limit L -> oo. We also fix £ = SU{2) = Sf. 

The Skyrme model cannot be solved analytically either in flat space, or on S^, apart 
from the case if L — 1 and B — 1 as mentioned above. However, in flat space there 
are analytic ansatze which give good qualitative and quantitative agreement with exact 
solutions obtained numerically. By ansatz we mean a test function that minimizes the 
energy within a given class of test functions. The lower the energy the better we expect 
the ansatz to approximate the exact solution. In the following we will generalize the 
rational map ansatz [J], which has been very successful in flat space, to S"£. 

A rational map is a holomorphic function from S 2 — > S 2 . Treating each S 2 as a Riemann 
sphere, with complex coordinates z and R, respectively, a rational map can be written as 

R{z) = 44 (3-1) 

where p(z) and q(z) are polynomials in z which are assumed to have no common factors. 
It has been shown by Donaldson |IT| , and also by Jarvis |18[], that there is a one-to-one 
correspondence between rational maps and monopoles. In flat space it has been found that 
many solutions of the Skyrme equation with baryon number B look rather like monopoles 
with monopole number equal to B. Therefore, rational maps can be used to approximate 
Skyrmions. 

A point x on S 3 is labelled by polar coordinates (//, 9, (/)) such that 

x = (L sin /i sin 9 cos 0, L sin /i sin 6* sin (f>, L sin cos 0, L cos /i) (3.2) 

= (L sin fi n(0, 0) , L cos /i) (3.3) 

where /i, 9 G [0,7r], and (f) G [0, 2iz] and h(9,<p) is the unit vector on S 2 . The 3-sphere can 
be thought of as a collection of 2-spheres with varying radius equal to Lsin/i. With the 
stereographic projection z = tan |e*^ the S 2 can be identified with a Riemann sphere using 
a single complex coordinate z. Alternatively, we can express the unit vector n in terms of 
z: 

ft W = T^VT2 (2Re(z),2Im(z),l - \z\ 2 ) . (3.4) 

1 + \Zf 

Similarly, points in the target S 3 can be labelled by (/, R) where / is an angular variable 
analogous to //, and R is a complex coordinate. The rational map ansatz simply states 
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that / = /(/i) and R = R(z). This is only consistent if sin/(/x) vanishes where sin /x does, 
i.e. /(0) = A^7r and /(7r) = N 2 it. The integer iVy = Ni — N 2 is a topological invariant 
and cannot be changed by deforming /(/x) smoothly In order to have a good limit for 
L — > oo we fix /(/x) such that i\r 2 = and set iVi = iV/. In analogy to flat space we 
define a Skyrmion to have Nf > 0, whereas for an anti-Skyrmion Nf < 0. Since Skyrmions 
and anti- Skyrmions are related by reflection we will only consider Skyrmions from now on. 
Therefore, our complete boundary conditions are: 

/(0) = N f n where N f > 



/(tt) = 0. (3.5) 

In contrast to the flat case these boundary conditions do not follow from a regularity 
argument. They are an artifact of our ansatz and have to be handled with care (see JTS| 
for a discussion). 

If we now use the notation of equation ( |2.10 ) we can write the Skyrme field in the 
following way: 

U(/i,z,z) = cos/(/x) + z sin /(/x) n(i?(z)) ■ T (3.6) 



exp(z /(/x) n(R(z)) ■ r) 



(3.7) 



where h(R(z)) is as in equation (|3.4| ). Equation ( |3.7|) looks quite similar to the well known 
spherically symmetric hedgehog ansatz [|IJ. In fact, the hedgehog ansatz corresponds to 
the special case R(z) = z. 

Now, we can apply the formulae of the previous section. The rational map ansatz gives 
rise to the following eigenvalues \\ of the strain tensor 



Ai 



L 

L 



and Ao = As 



sin / 1 + \z\' 
Lsin /x 1 + \R\ 



dR 



dz 



(3.1 



The minus sign in the expression for Ai is a consequence of our definition of positive baryon 
number in (|3.5| ). One advantage of the rational map ansatz is the decoupling of the radial 
and the angular strains. Starting from formula (|2.6|) the baryon number B can be written 
as a product of two integrals, one over /x and one over z and z in the following way: 



B 



(1 + 


z 


2 ) 


(l + \R 


2 ) 



dR 



dz 



2i dzdz 



1 + z 



2\2 



N f N 



R- 



(3.9) 
(3.10) 



The integral over \i is the integer Nf . The second integral is the pull-back of the area form 
on the target sphere of the rational map R(z). Taking the normalization into account this 
is just the degree Nr of the rational map R(z). In fact, it can be shown that if p(z) has 
the degree n p and q(z) has the degree n q then Nr 
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The energy can now be obtained from formula ( |2.4j ): 

sin 2 / 



E 



■ £12 / £12 



1 + 


2' 


2 


dR 


1 + 


/? 


2 


dz 



(3.11) 



+ 



sin 4 / 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



2i dzdz 



(l + \z 



2\2 



L 3 sin 2 \i d\i. 



Similarly to the baryon density, the integration over z and z and the integration over /i 
factorizes. Therefore, formula ( |3.11| ) can be rewritten as 



E = 4tt / [ f 2 L sin 2 fi + 2N R (C- + L) sin 2 / + X Sm \ f I d/x 

L Lsin /i, 



(3.12) 



where is the degree of the rational map, and X is the following special function on the 
space of rational maps: 



X 



1 

An 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



2i dzdz 



1+* 



2A2 ' 



(3.13) 



To minimize i? for a given baryon number B = NfNn, one should first minimize X over 
the space of rational maps of degree A^. This calculation was performed in 0, and the 
result is displayed in table |l| in the appendix. Then the profile function /(/i) is found by 
solving the following Euler-Lagrange equation: 



(2N R . 2 . 2 \ 2N R 2 

l-^-sm / + sin /i I + 2/ sm / ucos/i+ -^-/ sm/cos/ 

2X sin 3 / cos / 



- 2A^ R sin / cos / - 



L 2 sin 2 /x 







(3.14) 



where X now takes the constant value in table [I]. We require f(n) to be a solution of ( |3.14|) 
non-singular in [0, 7r]. Equation (|3.14| ) has to be solved numerically. It has regular singular 
points at the end points, i.e. close to these points the solution has the following power 
law behaviour: /(/i) ~ Nfir — A±fi p± for /i « and /(/i) ~ B±(n — fi) p± for /i « tt. Here 
P-t = |(±a/1 + 8N R — 1) and A± and 5-t are arbitrary constants. The solution /(//) is 
regular if the exponent is equal to p + at both end points. 

Equation ( |3.12| ) has an important discrete symmetry. The transformation 



/M - N f n - /(vr - n) 



(3.15) 



transforms a solution of ( p,14| ) into a solution which is also compatible with the boundary 
conditions (|3.5|) . Geometrically, /i — > tt — jj, is a reflection at the plane through the equator 
in physical space, whereas / — > Nfir — / is a reflection in target space. This means that 
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a solution which is localized for example at the south pole \i = 7r is transformed to a 
solution which is localized at the north pole fx = 0. Therefore, for fixed B there are two 
degenerate solutions unless the transformation ( 3.1 5|) is a symmetry of f(p) in which case 
there is only one symmetric solution. The symmetry ( |3. 15|) does not have an analogue in 
flat space. Also note that the transformations f(fx) — > /(tt — fi) and f(ft) — > —f(fx) take a 
Skyrmion with baryon number B into an anti-Skyrmion with baryon number —B. 

When we derived the energy of the rational map ansatz we used equation (|2.4j ). This 
could be transformed into equation ( |2.7| ) by "completing the square" . Having now calcu- 
lated the integral over the rational map, we can complete the square in a different way and 
re-express ( |3.12|) as 



E = AnL 



f sin fx + 



v/Jsin 2 { 
L sin fi 



+ 2N R sm 2 f 



r 



+ 1 



dfx 



-8n(2N R + VT) [ /'sin 2 /d/i. 
Jo 



(3.16) 



The second integral in ( |3.16| ) can be evaluated using f(0) = Nj-k and f(n) = 0. Since the 
first integral in (|3.16|) is positive the second integral provides us with the energy bound 



E > 4it 2 (V1 + 2N r )N j 



(3.17) 



> 12ir 2 N R N f . (3.18) 

The last inequality is valid because J > N R which can easily be shown by applying 
the Schwartz inequality. Therefore, inequality ( p. 17 ) "improves" the Faddeev-Bogomolny 



bound ( |2.8j ) where Vol(S) = 27r 2 and B = N R Nf. This rational map bound ( |3.17| ) is valid, 
if the fields obey the rational map ansatz. Exact solutions satisfy neither the rational map 
ansatz nor necessarily the bound ( |3.17| ). We will discuss this bound further in section [5[ 



Before ending this section we review one further ansatz. The symmetry group of S^, 
which is 0(4), contains an 0(2) x 0(2) subgroup. Jackson et al used this symmetry to 
obtain doubly axially-symmetric ansatze for Skyrmions of various baryon numbers B [[[2] . 
In order to make best use of the symmetry it is convenient to parameterize the 3-sphere 
with a different set of angles (x, «, (3). With these coordinates a point on S~[ can be written 
as 

x = (L sin x cos a, L sin x sin a, L cos x cos j3, L cos x sin (3) (3.19) 

where x e [0> f ] an d ct,j3 G [0,27r]. Now we can write an 0(2) x 0(2) symmetric field 
configuration U G SU{2) as 

U = sing cos pa + ir 3 sin g sin pa + vr\ cosg cos q/3 + zr 2 cos gsin q/3 (3.20) 

where p and q are integers and g = g(x) is a shape function. In this ansatz the eigenvalues 
A,- of the strain tensor (12. 21) are 



q' , v sin q , q cos q 

Ai = y, A 2 = f-^ and X 3 = j y -. 3.21 

L Lsmx L> c °s X 
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Using equation ( |2.4| ) we obtain the energy 



E = 2^L f ( g " + ^ + l^)M^h (3.22) 
7 V sm x cos z x 

+ 2 4~ I ( + ^0 + PV f V 2 ° SV ) ^(2x)dx 

L J \ \ sm x cos X / sm x cos X 



and the baryon number 



B = M /2 J 'si„ 9C o 8S dx. (3.23) 



The shape function g(x) can be calculated numerically by minimizing the energy E subject 
to the boundary conditions g(0) = and g(%) = f for various p, q and L. Note that these 



boundary conditions for g(x) imply that B = pq. In JL^] the energy for a given baryon 
number B was also minimized with respect to the radius L. Some of the results are 
displayed in figure [| of section ||. It was shown that the solutions are only stable as long as 
L is smaller than a 5-dependent critical length L cr i t .. Moreover, for large baryon number 
B the minimal energy for the optimal radius L is larger than the energy of B well separated 
Skyrmions in M 3 . Therefore, in this situation the ansatz fails badly. However, for small 
baryon number B and small radius L the ansatz is quite successful. For B = 1 it agrees 
with the known 0(4) symmetric solution. For B = 2 it predicts an 0(2) x 0(2) symmetric 
solution with a very low energy. Given that the exact B = 2 solution in flat space possesses 
an 0(2) x Z2 symmetry this configuration is likely to be the exact solution on S* 3 . 

One particular feature of configurations of the form ( |3.20| ) is that the following order 
parameter 0\ vanishes: 

O x = (a) 2 + (tt) 2 = (3.24) 

where ( ) means the average over the physical S^. It also turns out that configurations 
with p = q attain their minimal energy at particularly small values of L. Moreover, for 
p = q the solution has a symmetry similar to ( J3.15|) . The transformation 

g(x) - \ - g(\ - x) (3.25) 



transforms solutions of the Euler-Lagrange equation for the energy (|3.22 ) into each other. 
We will need these properties in section ||. 



4 The Shape Function as a Quasi- Conformal Map 

We discuss next the shape function f(fx) of the rational map ansatz. Since solving equation 
( p.!4| ) numerically provides little insight we derive an analytic shape function. This ansatz 
approximates the numerical shape function fairly accurately and also confirms geometric 
ideas. 

In H Manton approximated the shape function of the B = 1 Skyrmion by a conformal 
map. With this ansatz he showed that a delocalized B = 1 Skyrmion on a 3-sphere is 
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unstable for L > y/2. In this section we will generalize this idea for higher baryon number 
B. Using the rational map ansatz we will derive an ansatz for the shape function which 
is conformal in an average sense and which we will call quasi-conformal. In the following 
section, we will show that this quasi-conformal ansatz is a good approximation to the 
numerically calculated shape function. 

A Skyrmion is a map from physical space S\ labelled by (/i, z) to a target S\ = SU (2) 
labelled by (f,R). A map between two 3-spheres is conformal if their metrics only differ 
by a conformal factor: 

2i dz dz „, 2 / 2 . 2 „2i dR dR 



L dfi + L s ^, JT - W2= ^^Y{if' + ^f JTTm y (4.1) 

We are interested in fields which obey the rational map ansatz, i.e. f = f(fi) and R = R(z). 
Therefore, we make the following approximations to equation ( |4.1| ). Firstly, since /(//) is 
a function of /i only, we consider a conformal factor Q which is also only a function of 
/i. Secondly, we recall that according to equation ( |3.9| ) the integral over the target S 2 
is just An times the degree Nr of the rational map. Therefore, we replace the metric 
on the target S 2 by Nr times the metric on the physical S 2 . Locally, this is not a very 
good approximation, but averaged over the whole S 2 this reproduces the correct result. 
Since f(fi) is independent of z and z it can only detect the averaged value. With these 
approximations equation (|4.1|) can be simplified: 



r , , , t9 9 2i dz dz ^ . , 9 / , n9 9 2% dz dz \ 
L 2 d/x 2 + L sin 2 fi — = Q (jjt 2 df + sin 2 / N R . 4.2 

(1 + zz) 2 \ (l + ZZ) 2 J 

Eliminating O(yu) from equation ( f4.2[ ) we obtain the following differential equation for /(//): 

(f(n)) 2 sm 2 v = N R sm 2 f(ri. (4.3) 

We are interested in solutions which obey the boundary conditions ( |3.5|) . It is convenient 
to replace sin/(/i) by sin(7r — f(fj)) before taking the square root: 

/V) sin/i = ±v / A^sin( 7 r - (4.4) 



For the negative sign the solution of equation (O) diverges at the boundary and can 



therefore be discarded. Equation ( |4.4j ) is solved by separation of variables, and we obtain 



/(/f) = 7r — 2 arctan ^tan j 



(4.5) 



where k is a positive integration constant. A negative k would lead to a negative baryon 
number and is incompatible with the boundary conditions (|3.5| ) whereas k = is just the 
trivial solution with baryon number B = 0. The quasi-conformal shape function (|4.5| ) 
satisfies our boundary conditions ( |3.5| ) if and only if A^y = 1. With equation ( |3.10| ) we 
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obtain B = Nr. The shape function (|4.5|) has the following important property. Using the 
identity 

1 7T 

arctanx + arctan — = — for x > (4-6) 
x 2 



it is easy to show that for k = 1 the shape function ( |4.5| ) is invariant under the discrete 
symmetry (|3.15| ): 

/(Ai)=7r-/(ir-/i). (4-7) 

This means that for k = 1 the solution is neither localized at the south pole nor at the north 
pole. For Nr = 1 the solution is delocalized over the whole sphere whereas for Nr > 1 it 
is rather localized at the equator fj, — |. By Taylor expanding equation ( |4.5| ) near /i = 
we find that the shape function /(/i) behaves like /(/i) ~ vr — 2k(^)^R. Recall that for 
the exact solution of the Euler-Lagrange equation (|3.14| ) we obtained f(fi) ~ n — A + /j, p+ 
where p + — — | + \\/^ + 8N R . So at n = the shape function ( fl5l) and the exact solution 
have a similar behaviour. The same is true for the other boundary = ty. 

Now, we can substitute the shape function ( [4.5[ ) into the equation for the energy ( |3.12| ) 
and obtain 

E = 16tt (?>LN R h{k) + | (2N R 2 + X) I 2 {k)j (4.8) 
where 

Jo (1 + fc 2 y 2V I 1 + y 2 ) 

J 2 (fc) = / ^ S^dy. (4.10) 

1 ; io 2(1 + A;V^) 4 

In (f4.9|) and ( |4.10| ) we have simplified the integrand by substituting y = tan|. The 
integrals h{k) and hik) can be evaluated in closed form if \/Nr is an integer. Therefore, 
we will concentrate on Nr = 1,4 and 9 in the remainder of this section. 



and 
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Figure 1: The quasi- conformal shape function f(fi), its derivative f'(^), and the averaged 
energy density E(fi) for B = 1 and 5 = 4. The parameter /c takes the values 4, 1,4. 



Figure [I] shows the shape function ( [4.5| ) and its derivative for S = iV^ = 1,4 and for 
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k = |, 1, 4. Figure [I] also shows energy density E(fi) averaged over the 2-sphere, which is 
the integrand of G3.12|) . The shape function (^4.5|) has symmetry ([4.7|) for k = 1. For k > 1 
the shape function is localized around the north pole /i = 0, whereas for k < 1 it is localized 
around the south pole \i = n. This is particularly obvious for its derivative f\n). Figure 
[I] also illustrates that for higher baryon number B the symmetric solution becomes more 
and more localized at the equator \i — |. Furthermore, the energy E(/i) has been plotted 
in order to compare the result to the usual Skyrme model. E(fj,) corresponds to an average 
radial density. When the energy is peaked close to fx = this is the usual Skyrmion. If the 
energy is centered around \i = | this corresponds to a shell-like configuration. Solutions 
centered around /i = ir will correspond in the limit L — > oo to a configuration centered 
around infinity, with an infinite energy. In order to calculate E(n) in figure |1] we have 
to know the values for X and L. As we shall see in the following, for a given B, k ^ 1 



minimizes the energy ( |4.8| ) for a unique radius L. k = 1 is a solution for many different L, 
but among these solutions L crit . takes a special place. For B = 1, E(fi) has been calculated 
by setting X = 1, and L = 2.21 corresponding to k = 4 and fc = t. For = 1 we 
have displayed for the critical length L CT u. = 1.41. Similarly, for B = 4, L = 2.47 

corresponds to k = 4 and A; = |, and we have also displayed for L criti = 2.09. X takes 
the value X = 20.65 of table |T] in the appendix. 

We will now discuss the energy (|4.8| ) as a function of L. For 5 = 1, and hence Nr = 1, 
the rational map is the identity map R{z) = z. The integral ( |3.13| ) can be evaluated 
explicitly. We obtain X = 1 so that the energy ( |4.8| ) is 



E = 127T 2 ( 2 f + — (k + -)). (4.11) 
\k + l + 2 AL\ kJJ y } 

It is convenient to set a = k + | and to calculate the minimum of the energy with respect 
to a. E has a single minimum at a = VSL - 2. For L < \/2, 

cto is not attainable for real 

fc, so that the minimum energy occurs at a = 2 and k — 1. In fact, up to a minus sign this 
is the identity map between the physical and the target Sf. The energy E can now be 
written as 

£ = 6tt 2 (l + ^. (4.12) 
For L > \[2 the minimum occurs where 

k + l = V8L-2. (4.13) 
Equation (4.13) has two solutions. They are related by the symmetry ( |3.15| ) and correspond 

to solutions loraliwH at onp of the nolfis Thfiir prmrcv is 



to solutions localized at one of the poles. Their energy is 

E=127r 2 (V2-^j (4.14) 

which is lower than ( |4.12|) . In the limit L — > oo the energy becomes E = 127t 2 a/2, and the 
shape function turns into f(r) = 2 arctan(v^8^) where r = L/i. 
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Figure ^| shows the energy E of the Skyrmion as a function of L. For small L the 
energy scales like 4. At the optimal value of the energy L opt , = 1 the energy reaches its 
minimum E opt . = 12n 2 . The Skyrmion is symmetric if the radius L is below the critical 
radius L cr i t . = \[2. For L > L cr i t . there are two Skyrmions with the same energy. This 
phenomenon is often called bifurcation. There is still a symmetric solution but it is no 
longer stable. For L —>■ oo the energy tends to the K 3 value. For B > 1 the energy E 
shows similar behaviour. The values L opt , L crit and E opt depend on the baryon number 
B and are good numbers to characterize the behaviour of the solutions. 




0.5 1 1.5 2 2.5 3 3.5 4 

Radius L 



Figure 2: The energy E in units of 12ir 2 of a B = 1 Skyrmion as a function of the radius 
L. 



Now, we consider Nr = 4. In this case the optimal rational map has octahedral 
symmetry and can be written in the following form: 

— 2y 6iz z + 1 

We obtain I w 20.65. Using expression ( |4.8|) the energy is 

_ 24n 2 V2L(f3 3 -6(3 + 4V2) 5^(32 + 1)(3 

E = W^W + lei (416) 

where (3 = \fk + i plays the same role as a for S = 1. We can now minimize the energy 
with respect to (5. This gives rise to the cubic equation 

(P - V2) 3 - P(J3 - \/2) - Q = 0, (4.17) 
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where P = ^^i) anc ^ Q = ^§2+ff • There is a real solution (3q for all L but its exact form 
is long and not very illuminating. For L < L CT u. we have (3$ < 2 so that j3 — 2 and k — 1 is 
again the minimum, and the solution is localized around the equator. The critical length 



can be calculated by substituting (3 = 2 into equation (4.17), and we obtain 



1 / 15(32 + J) 



■'crit. 



24 V 8V2- 11 

For k — 1 the energy simplifies: 



2.091 (4.18) 



g = 24^(2-y/2) + 5 ^ 2 <f + I) . (4.19) 
The optimal length can be calculated by minimizing E with respect to L: 



For L > L crit there are two solution such that yk + -4= = /3 related by the symmetry 
( |3.15| ). The corresponding expression for the energy can be derived explicitly by setting 
(3 = j3 in equation ( |4.16| ) but it is rather lengthy. 



The limit L — > 00 can be derived from equation ( |4.17 ). We are interested in the solution 



where (3 is of order L such that the energy E in equation ( |4. 16| ) is finite. Therefore, we 



can neglect the last term in ( f4.17| ), and the a/2, and obtain (3 = yP. Inserting (3 into the 



expression for the energy (|4.16|) we obtain = 4.684 which corresponds to 1.171 per 
baryon. 

For Nr = 9 the optimal rational map has D^d symmetry and X pa 109.3. It is again 
possible to introduce an auxiliary variable 7 = ks + k~3 . However, minimizing the energy 
leads to a polynomial of degree 5 in 7. Therefore it is no longer possible to calculate the 
solution explicitly. But the general behaviour remains the same. The solution is symmetric 
below a critical length L crit . which is determined by substituting 7 = 2 into the equation 
for minimizing the energy. We obtain 



L cnt = V323190 + 1195 J pa 2.887. (4.21) 

513 



The optimal radius can also be calculated analytically: 

1 



L op t. = — V374220 + 2310 J pa 2.683. (4.22) 
297 

For L > L crit . there are again two degenerate solutions. In the limit L —>■ 00 we can balance 
the terms which are of highest order in L. This leads to 7 ~ 0.779L, and the normalized 
energy is = 10.274 which is 1.142 per baryon. 
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In the last part of this section we will discuss the limit L — > oo for general Nr. Setting 
r = L/j, the quasi-conformal ansatz (|4.5|) tends to 




/(r)=7r-2axctanN— j J (1.23) 
where i?o is a free parameter. The energy can be written as 

E = t (^f Jl + (2iV * + J) Roh ) (424) 

However, now the integrals I± and I2 are independent of Rq and depend only on Nr: 

poo r 2 v / 7Vfl 

h = / ^^dr (4.25) 

Jo (l + r 2 v^)2 1 > 

h = / ^^dr (4.26) 

It is straightforward to minimize the energy (|4.24|) with respect to Rq. In figure ^| and ^, 
in the following section we have evaluated the minimal energy E(Rq = i2o,mm)> an d the 
optimal radius Ro, m in, respectively. The parameter -Ro.mm has a natural interpretation as 
the size of the Skyrmion: For r = -Ro.min the shape function f(r) has reached the value |. 
Furthermore, when the baryon density is integrated over a ball of radius i?o,mm, then this 
is just half the total baryon number: 

/ sin f df = —. (4.27) 

7T Jo 2 

It is also worth noting that equations ( |4.5|) for k = 1 and ( |4.23| ) are related by r = 
-Ro,mmtan|. This is a stereographic projection from a 3-sphere of radius Ro,min to R 3 . 
Minimizing the energy we see that a Skyrmion in flat space is related to a Skyrmion on 
the 3-sphere with an optimal radius L opt . 



5 Numerical Results 



In this section we discuss the numerical solution of equation (|3.14 ) for the shape function 



f(fji) of the rational map ansatz and the properties of the resulting fields, and we compare it 
to the analytical solutions from the previous section. Here, we only consider configurations 
with Nf = 1 such that B = Nr. A configuration with Nf = 2 is discussed in the following 
section. 

The numerical solution of ( |3.14| ) is calculated using a relaxation method. For symmetric 
initial conditions we obtain the symmetric solution which is a saddle point for L > L cr i t .. 
For asymmetric initial conditions we obtain the minimum energy configuration, i.e. the 
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Skyrmion. In figure [5] we compare the numerical shape function with the quasi- conformal 



for B = 4 and B = 9 at the optimal radius L = L opt . The optimal radius L. 



J opt. 



ansatz 

implies k — 1 in the quasi- conformal ansatz. The numerical result and the quasi-conformal 
ansatz show good agreement, in particular in the region around the equator jj, — |. This 
region is most important for the energy ( p,12| ) since both (/'(aO) 2 an d sin 2 /(/x) are large 
for fi fa |. Near the poles the numerical solution is less steep than the quasi-conformal 
ansatz as we expect from their exponents: p + > \f~B. 




Figure 3: Comparison of the numerical shape function with the quasi-conformal ansatz for 
B = 4 and B = 9. 



In figure § we show the energy E as a function of L for B — 1 , . . . , 4 and compare 
it to our analytic calculations of the previous section. In all the figures the energy E 
is normalized by 12ix 2 so that the Faddeev-Bogomolny bound ( |2.8|) is equal to \B\. For 
B = 1 we plot energy ( [4.12| ) of the symmetric solution which gives the exact result. For 
L > L crit . = \/2 we also plot the energy (|4.14|) . Our ansatz predicts the correct critical 
radius but energy ( |4.14| ) is slightly too high. For B = 2 and B = 3 we only display the 
energies of the symmetric solutions using formula (|4.8|) and setting k = 1. The integrals 



/i(l) and i2(l) in ( }4-8|) are calculated numerically. Again there is very good agreement for 
small L, and for L > L cr i t . the ansatz is also very close to the symmetric solution. For B = 4 
we plot the analytic expression for the symmetric solution (|4.19 ). For L > L q r '~^ con ^' w 2.09 
we also plot (|4.16|) with (3 = 0o 



The critical radius L q c ~^ on ^' is slightly higher than 
Lcrit enc = 2.071, and the energy is also slightly too high. In figure f| we also mark the 
minimal energy E™£ 

.0(2) 



numeric 



for the optimal radius L 



numeric 
opt. 



and for comparison the optimal 



energy E Q p t for the doubly axially-symmetric ansatz at its optimal radius L 



S 2) - We will 
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discuss these quantities in figure [IJ 

In summary, all the analytical results are in good agreement with the numerical solution 
for a large range of L. This is quite remarkable given that the ansatz for the shape function 
only depends on one parameter k. 




In figure [5] we compare the optimal energy of the quasi-conformal map E q opt con with the 

energy of the doubly axially-symmetric ansatz £w for all baryon numbers B = 1, . . . , 9. 
For B = 1 both reproduce the exact result. For B = 2 the doubly axially-symmetric ansatz 
is believed to be the exact solution because its symmetry is compatible with the results in 
flat space. Yet, for B > 2 the energy E q ~ t con ^' is lower than E° p f\ 

Figure [5] also shows that the quasi-conformal ansatz always has a lower energy than the 
numerical solution in M 3 . We conclude that L opt . for the exact solution will be finite. The 
horizontal line at E = 1.232 in figure corresponds to B well separated single Skyrmions. 
For B > 8 the energy of the doubly axially-symmetric ansatz is above this line which is 
unphysical. Note that the energy per baryon in flat space Eyp decreases as we increase the 
baryon number. We would obtain a similar behaviour if we kept the radius L fixed. Yet, 
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Figure 5: The energies E^ con ^\ E^\ E®3 and EBound per baryon as a function of the 
baryon number B. 



we then have to make the radius L large enough so that the largest Skyrmion "fits". The 
physical interpretation is that one multi-Skyrmion is more stable than a number of smaller 
Skyrmions. 

We also display the rational map bound E Bound ( |3.17|) . For B = 2 the energy of the 
doubly axially-symmetric ansatz is lower than EBound- This confirms our warning that 
energies of exact solutions do not have to be greater than EBound- However, figure [| shows 
that the bound follows the general behaviour of E q ~ t con ^' and also of E^z. 

In figure ^| we compare the energy of the quasi- conformal map E^°° in the limit 
L — > oo to the rational map energies E™™ enc and the energy E-^3 of the exact solution 
in flat space. E^ con ^' is calculated by minimizing the energy of expression ( [4.24 ) with 
respect to Rq. For B = 1,4,9 the energy E q ~ t ^ on ^' agrees with the results of section [|. 
Figure ^| shows that the quasi-conformal ansatz approximates the numerical solution well. 
The relative difference between E'^ t " m and £^ me " c is monotonically decreasing and is 
less than 3% for B = 4. This means that the error of the rational map ansatz with a 
numerical shape function and the quasi-conformal map ansatz are of the same order of 
magnitude. 

In figure [7] we compare the different optimal radii. The optimal radius of the quasi- 
conformal ansatz L q ^ con is always slightly smaller than L"p t mer * c . In all cases the value of 
the critical radius L crit _ is greater than the optimal radius so that the results are consis- 
tent. The critical radius can be calculated in two different ways. In the first approach we 
determine the point at which the solution bifurcates. We use the fact that the symmet- 
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as a function of B. 
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ric initial condition relaxes into the symmetric solution which is always present, whereas 
an asymmetric initial condition only relaxes to the symmetric solution when there is no 
asymmetric solution. This method leads to rather large error bars because it is difficult to 
decide whether two numerically calculated solutions are different. The second approach is 
to consider the value of the following integral: 

(/M+Zfr-^-Tlfd/X. (5.1) 



This integral is zero if the solution possesses the symmetry ( 3.15|) and is nonzero other- 
wise. We determine the value of L where the integral ( |5.1| ) first becomes nonzero. Both 
approaches give the same results, but the second one is much more accurate. At the crit- 
ical value there is a jump in the value of the integral (|5.1|) . This is also a verification of 
our analytic result that the solution is symmetric for L < L cr i t .. As mentioned before, for 
B = 1 the critical radius is in agreement with the analytic result L cr i t . = a/2. For B = 4 
the critical value of the quasi-conformal ansatz L q c ~ i ^ on ^' = 2.091 in ( 4.18Q is only slightly 
higher than the critical radius of the rational map L™™ ertc = 2.071. Similarly, for B = 9 
we calculated L^ con/ - = 2.887 in flPip whereas = 2.866. 

In figure [7| we also find that Ro,min is always smaller than L q ~ t con ^' . The difference 
between these lengths becomes smaller the larger the baryon number is. This confirms our 
interpretation stated at the end of the previous section, that the Skyrmion in flat space is 
related to the Skyrmion on a 3-sphere with optimal radius. 



6 Phase Transitions 

In this section we will describe phase transitions on S 3 as L is varied. There is no numerical 
work solving the full set of coupled partial differential equations for £>£ as there is in flat 
space. Therefore, we have to rely on our ansatze. First, we will discuss two different order 
parameters. Then we will describe particular situations in more detail. 

In R 3 the average (a) is equal to 1 because of the boundary conditions at infinity. 
However, on it is easy to show that (a) vanishes for L < L cr u. when the solution is 
symmetric with respect to south and north pole. Therefore, we can choose 0\ = (a) 2 as 
an order parameter.^ It has the necessary property that it is nonzero above the phase 
transition L > L cr i t . and identically zero below the transition. As we already pointed out 
in section |3| this order parameter also vanishes identically for all doubly axially-symmetric 
configurations. This is consistent with the fact that for L > L cr i t . these configurations are 
bad approximations. Another possible order parameter would be O2 = (cr 2 ) — 7. This 
parameter O2 is motivated by the fact that enhanced chiral symmetry, that is a symmetry 
that mixes a and 7Tj, could result in a vanishing of C>2- 

For B = 1 the situation is well understood. Both parameters 0\ and O2 vanish at the 
same time. Physically, this means that the localization-delocalization transition occurs at 

2 A more symmetric choice is (ct) 2 + (vr^) 2 , Q. 
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the same time as the restoration of chiral symmetry ||. In fact, for L < y/2 the solution 
possesses full chiral symmetry 50(4). For B > 1 the chiral symmetry will be a subgroup 
of 50(4). 

Next we discuss B = 4 and use both the rational map ansatz and the doubly axially- 



symmetric ansatz. In figure |8(a)| we plot their energies. For L > L crit . the Skyrmion is 
localized at one of the poles. For smaller L the Skyrmion is localized around the equator. 
This is the first phase transition and the order parameter is 0\. However, there is a 
radius L cr i t / where the doubly axially-symmetric ansatz with p = 2 and q = 2 becomes 
the minimal energy solution. As pointed out in section |3| if p — q the solutions have 



the symmetry ( |3.25| ). For this special case we checked numerically that the symmetry is 



actually a symmetry of the shape function g(x)'- 

#(x) = |-#(|-x) (6- 1 ) 

This gives rise to a symmetry between a and tt\ and leads to the vanishing of the order 
parameter O2. Therefore, we have a second phase transition at L = L crit .i. 




(a) The energy of B = 4 Skyrmions with (b) The energy of two B = 8 Skyrmions, one 

the rational map ansatz and with the doubly with Nf — 1 and one with Nf = 2. 

axially-symmetric ansatz. 



Figure 8: Phase transitions in the Skyrme model 



As a last example we describe a Skyrmion with baryon number B = 8. For large L 
the Nf = 1 solution has minimal energy. However, there is a saddle point solution where 
two B = 4 Skyrmions are located at opposite poles of the 3-sphere. If L is sufficiently 
small there is a solution with Nf = 2 which has lower energy than the one with Nf = 1 as 
shown in figure |S(b)| .P| This configuration is particularly interesting because both 5 = 4 

3 Table || in the appendix gives an indication when to expect a configuration with Nf > 1 to be more 
stable than the corresponding configuration with Nf = 1. 
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Skyrmions have cubic symmetry. For sufficiently small radius Lh c these two Skyrmions 
could combine and form a solution with a hypercubic symmetry. This symmetry would be 
strong enough to restore the symmetry between the a- and 7Tj-fields. The order parameter 
O2 would vanish. In this scenario there would be a localization-delocalization transition 
at L cr it., whereas the chiral symmetry is further enhanced at Lhc- Unfortunately, these 
discrete subgroups of the full 5*0(4) symmetry cannot be studied with the rational map 
ansatz. 

In physical terms we propose the following picture. Reducing the radius L increases 
the baryon density. If we start with a B > 1 Skyrmion at large length this is well described 
by the rational map ansatz. There is a phase transition at the critical radius L cr u. where 
the Skyrmion becomes delocalized over the 3-sphere and the order parameter 0\ vanishes. 
Then, one or more phase transitions follow which make the configuration more and more 
symmetric, and chiral symmetry is partially restored. Note however, the Skyrme model is 
an effective theory which is only valid for low energies. If the radius is too small we expect 
that further terms in the effective action become physically important. 



7 Conclusion 

In this paper we have described the rational map ansatz applied to Skyrmions on a 3-sphere 
of radius L. We have calculated the energy E of a Skyrmion as a function of L for baryon 
number B — 1, ... ,9 and found the following behaviour: For small L the energy E scales 
like 4-. E has a global minimum at an optimal radius L opt .. When the radius L is further 
increased there is a bifurcation point at a critical length L crit : below L crit the Skyrmion 
is symmetric under reflection at the plane through the equator of S 3 . Above L crit , there 
are two degenerate Skyrmions one of which is localized at the north pole whereas the other 
one is localized at the south pole. For L — > 00 the ansatz tends to the rational map ansatz 
for M 3 , H]. For B = 1 and L = 1 our ansatz reproduces the known exact solution \TB 



For B = 2 our results are worse for small L than the doubly axially-symmetric ansatz in 
|T2|]. However, for B > 2 the energies of our solutions are lower than any solutions known 
to date, including the doubly axially-symmetric solutions. 

We have also derived an analytic expression for the shape function by imposing that 
the metrics on the physical S 3 and on the target S 3 are conformal in an average sense. 
This quasi-conformal ansatz for the shape function depends only on one parameter k. By 
varying k we obtain very good agreement with the numerical results for a large range of L, 
even for flat space. It is worth noting that the relative error of the quasi-conformal ansatz 
decreases with increasing baryon number. We have shown that the solution is localized at 
the equator if k = 1. For B = 1, 4, 9 we have shown that the Skyrmion really is symmetric 
below a critical value L cr u, i.e. k = 1 gives the minimal energy solution. 

One particular property of the quasi-conformal shape function is that it becomes more 
and more localized around the equator as the baryon number increases. To be more precise, 
the derivative of the shape function, which is connected to the baryon density, has a peak 
at the equator. Since this ansatz was derived with the assumption that the map between 
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the metrics is "as conformal as possible" we have found a geometric explanation of why 
Skyrmions on the 3-sphere are shell-like. In flat space the map cannot be conformal. 
However, one "half" of the Skyrmion in flat space resembles "half" a Skyrmion on a 
3-sphere with optimal radius. This fits well with the observation in section |5| that the 
parameter Ro,min which is a measure of the size of the Skyrmion agrees well with the 
optimal length /,™« i me " c . 

This line of thought can be carried even further. Skyrmions on the 3-sphere might be 
a reasonable model for nuclei, once they are quantized. The main advantage of this model 
is that on the 3-sphere one-loop corrections are expected to play a far less important role 
than in flat space. In particular, if the Skyrmion lives on the 3-sphere with the optimal 
radius we expect these corrections to be so small that the predictions of the Skyrme model 
could be compared with experiment. The main motivation for this claim are the promising 



results in |2(J for the B = 1 Skyrmion on the 3-sphere, and the fact that they predict a 
different value for the Skyrme coupling which is in agreement with one-loop calculations 
performed in M and the Skyrme coupling therein. 
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Appendix 

A.l 1 for N = 1,..., 9 and 17 



Table [I] shows the minimal value of the integral X in equation ( |3.13| ) as a function of Nr, 



as has been calculated in |7j . "APPROX" is the energy of the rational map ansatz for flat 
space Skyrmions where the corresponding variational equation for the shape function has 
been solved numerically. The column "TRUE" shows the value of the energy obtained by a 
numerical solution of the full set of partial differential equations. The next column "SYM" 
shows the symmetries of the Skyrmions. The numerical solutions and the rational map 
ansatz give the same symmetry. Note that the respective values for B = 9 in table [l] are 
taken from [2T|. There is evidence that for B = 9 the symmetry and the other respective 



values in the table are different from the ones cited in 0. 

The last two columns are the optimal length L opt . and the optimal energy E™mencai 
which are displayed in figure [7| and || respectively. As described in chapter ||, E™memcai 
has been obtained by using the rational map of table [l| and calculating the shape function 
numerically. 

Following |7j we also include the values for B = 17. In this case the solution in flat 
space has marginally lower energy than £;™« me ™ c . However, we still expect that the rational 
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map ansatz is a good approximation even in this situation. 



Nr 


1 


APPROX 


TRUE 


SYM 


Lopt. 


j^numeric 
opt 


1 


1.00 


1.232 


1.232 


0(3) 


1.000 


1.000 


2 


5.81 


1.208 


1.171 


0(2) x Z 2 


1.393 


1.078 


3 


13.58 


1.184 


1.143 


T d 


1.669 


1.093 


4 


20.65 


1.137 


1.116 


o h 


1.829 


1.065 


5 


35.75 


1.147 


1.116 


D 2d 


2.074 


1.089 


6 


50.76 


1.137 


1.109 


Dm 


2.250 


1.088 


7 


60.87 


1.107 


1.099 


Y h 


2.346 


1.064 


8 


85.63 


1.118 


1.100 


D 6d 


2.544 


1.080 


9 


109.3 


1.098 


1.093 


D M 


2.696 


1.083 


17 


367.41 


1.092 


1.073 


Yh 


3.614 


1.076 



Table 1: Values of X, the energies in flat space (APPROX, TRUE), the symmetry (SYM), 
and the optimal length L opt . and energy £ , ™« merjc Q n S 3 . All values are given for baryon 
number B = Nr for B = 1, . . . , 9 and B = 17. 



A. 2 The Rational Map Bound for Different Splittings B = N f N R 

In table we compare the rational map bound for different Nf. As the rational map bound 
is just the topological bound for Nr = 1, we have omitted Nr = 1 from table |2|. 



(Nf, Nr) 


(1,4) 


(2,2) 


(1,6) 


(2,3) 


(3,2) 


(1,8) 


(2,4) 


(4,2) 


(1,9) 


(3,3) 


D Bound 


4.181 


4.274 


6.375 


6.457 


6.410 


8.418 


8.363 


8.547 


9.541 


9.685 



Table 2: The value of the generalized Faddeev-Bogomolny bound (|3.17| ) for various values 
of N f and N R . 

Table |2| can be used to estimate which combination of Nf and Nr might be more stable 
for a given baryon number B. It turns out that for B = 8 the bound is lower for Nf = 2 
than for Nf = 1. It is also probable that the (Nf = 3, Nr = 2) Skyrmion has lower energy 
than the (Nf = 2, Nr = 3) Skyrmion. 
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